OTIC  FILE  COPY 


AFOSR.TR.  8  6-0436 


OPERATIONS  RESEARCH  AND  SYSTEMS  ANALYSIS 


in 

h* 

CM 


< 

I 

Q 

< 


Markov  Algorithms  for  Computing  the 
Reliability  of  Staged  Networks 

R .  C .  H .  C heng 

Technical  Report  No.  UNC/ORSA/TR-86/8 
April  1986 


UNIVERSITY  OF  NORTH  CAROLINA 
AT  CHAPEL  HILL 


DT1C 

ELECTFPf* 
JUL  2  8 1986  I 


86  7  23  226 


Approved  for  public  release  • 
distribution  uni  United. 


Markov  Algorithms  for  Computing  the 
Reliability  of  Staged  Networks 

R  .  C  .  H .  Cheng 

Technical  Report  No.  UNC/ORSA/TR-86/8 
April  1986 


Curriculum  in  Operations  Research 
and  Systems  Analysis 

University  of  North  Carolina  at  Chapel  Hill 


l 

s 


DT1C 

ELECTEj 
JUL2  8 1986 


D 


This  research  was  supported  by  the  Air  Force  Office  of 
Scientific  Research  under  grant  A FOSR - 8 4- 01  4 0  .  Reproduction  in 
whole  or  part  is  permitted  for  any  purpose  of  the  United  States 
government . 


distribution  statement  n 

Approved  lor  public  releoMl  I 
Piitttbuhoo  Unlimited 


ATR?”-  r  rri'TFTC  RESEARCH  (AFSC, 

1  C"  :?/N7V!TTSL  TO  DTIC 
Thin  technical  report  h?.n  been  reviewed  and  is 
approved  for  public  "eleass  IAW  AFR  190-12. 
Distribution  in  unlimited. 

MATTHEW  J.  KEEPER 

Chief,  Technical  Information  Division 


t 


ABSTRACT 


Certain  commonly  occurring  types  of  network,  whether 
directed  or  undirected,  exhibit  a  staged  structure.  Two 
algorithms,  based  on  node  partitioning,  are  presented  which  take 
advantage  of  such  structure  and  which  use  a  Markov  transition- 
probability  form  of  recursion.  The  algorithm  for  directed 
networks  is  related  to  the  Markov  chain  formulation  of  Bailey  and 
Kulkarni,  but  for  undirected  networks  a  more  detailed  form  of 
state  definition  is  used  related  to  one  suggested  by  Rosenthal. 

The  computational  advantages  of  the  algorithms  are  discussed 
and  some  numerical  results  presented. 


Key  Words:  Network  Reliability,  Source-to-sink  connectedness, 
Node  Partition,  Staged  Network,  Recursive  Algorithm. 


1  .  INTRODUCTION 


Consider  a  network  G  with  node  set  N  and  arc  set  A.  The 
nodes  are  perfect  but  arcs  fail  randomly  and  independently  with 
probability  1-pa,  aeA.  Let  s  and  t  be  two  specified  nodes.  The 
calculation  of 

Pst  *  prob  (there  is  a  path  from  s  to  t)  (1.1) 
is  called  the  s-t  connectedness  problem. 

This  problem  is  NP  hard  [10].  However,  many  algorithms 
exist  which  solve  it  by  enumerating  key  structures  in  the  network. 
If  the  number  of  these  is  small  then  such  an  algorithm  may  be 
computationally  efficient.  Examples  include  spanning  trees  [2], 
acyclic  subgraphs  [8],  K-graphs  [7]  and  cutsets  [ U  ] . 

A  number  of  methods  use  node  partition  [1,6,9].  Though  none 
of  these  make  any  special  assumption  about  the  form  of  the  net¬ 
work,  node  partition  seems  to  be  a  particularly  good  approach  for 
networks  which  have  a  staged  structure.  In  this  paper  we  consider 
two  node  partition  algorithms  which  exploit  such  structure. 

Our  starting  point  is  a  decomposition  used  by  Shogan  [9]  for 
directed  networks  where  the  nodes  are  divided  into  disjoint 
groups  -  which  we  call  stages  -  and  ’events'  are  then  defined  at 
each  stage  in  a  way  which  allows  their  probabilities  of 
occurrence  to  be  calculated  from  those  of  the  previous  stage. 
Shogan  describes  an  algorithm  based  on  path  enumeration  and 
mentions,  but  does  not  follow  the  possibility  of  one  based  on 
state  enumeration.  We  suggest  a  modification  of  the  definition 
of  event  given  by  Shogan  which  gives  an  algorithm  of  this  latter 
type  having  a  simpler  form  than  that  based  on  path  enumeration. 


Undirected  networks  can  often  also  be  decomposed  into  stages 


Our  main  result  will  be  to  show  how  the  algorithm  can  be  adapted 
to  deal  with  this  case,  too. 

Both  versions  of  the  algorithm  are  related  to  ones  already 
cited.  The  directed  network  version  is  related  to  that  proposed 
by  Bailey  and  Kulkarni  [1],  whilst  the  undirected  network 
algorithm  is  a  recursive  version  of  the  general,  but  not  fully 
specified,  framework  described  by  Rosenthal  [6]. 

As  well  as  describing  the  algorithms,  we  assess  their 
computational  advantages  -  and  drawbacks  -  including  some 
numerical  results. 

For  clarity,  the  bulk  of  our  discussion  assumes  a  form  of 
stage  decomposition  that  is  simpler  than  necessary.  In  section  5 
we  show  how  the  stage  decomposition  can  be  generalized  to  allow 
versions  of  the  algorithms  that  are,  in  certain  situations,  much 
more  efficient  than  the  prototype  versions. 

2.  DEFINITIONS  AND  NOTATION 
2.1  Staged  Networks 

Let  G  be  a  network  with  node  set  N  and  arc  set  A.  An  arc 
connecting  node  i  to  node  j  is  denoted  by  (i,j).  We  shall  say 
that  G  is  a  staged  network  if  N  can  be  divided  into  a  number,  Q 
say,  of  disjoint  sets 

S1  -  i s } .  S2 , . . . , Sq  =  { t  }  (2.1) 

such  that 


-3- 


S:  if  (i,j)eA  and  jeSq  then  ieSq_iuSq.  (2.2) 

(Shogan  gives  a  slightly  more  general  version  which  allows 
ieSiuSq_iuSq  in  (2.2).  The  generalized  definition  of  Section  5 
incorporates  this  case.) 

Condition  S  ensures  that  any  path  from  s  to  t  can  only  go 
from  nodes  of  one  stage  to  those  of  the  next  stage.  Thus  no 
'backtracking'  from  a  higher  to  lower  stage  is  possible. 

The  node-stages  induce  a  corresponding  stage  structure  on 
the  arcs  A,  dividing  these  into  disjoint  sets: 

A2*  A  ^ . Aq 

where 

A  q  *  1 (  i  » j  )  |  ieSq_iuSq  and  jeSq}.  (2.3) 

To  note  which  arcs  of  Aq  are  up  (operating)  and  which  are 
down  (failed),  we  define 

Xj  -  1  or  0  according  as  arc  ieAq  is  up  or  down,  and 

call 

x  -  (X1  ,X2 . X|  A  |  ).  (2.i») 

a  failure  pattern  of  Aq.  It  will  be  clear,  from  the  context, 
which  subset  Aq  is  being  referred  to,  so  the  dependence  of  X  on  q 
will  be  suppressed. 

Those  arcs  which,  under  X,  are  up,  will  be  called  the  set  of 
up-arcs  of  Aq.  The  probability  of  obtaining  X  is 

n  P  n  (1-P  >. 

a : X  - 1  a:X  -0  a 
a  a 


p(X) 


(2.5) 
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I'  Ql 

There  are  2  M  distinct  failure  patterns  of  A  and  we  denote 
the  set  of  all  these  as  Fq. 

Consider  now  an  undirected  network  G.  We  shall  say  that  G 
is  staged  if  the  nodes  can  be  divided  into  disjoint  sets 

Si  -  { s } ,  S2 ..... Sq  -  { t  } 

such  that 


T:  if  (i,j)  eA  and  ieSq,  1<q<Q  then  j eSq_ 1 u Squ Sq+ 1 .  (2.6) 

Condition  T  requires  that  an  arc  can  only  join  two  nodes  of 
the  same  or  adjacent  stages.  Decomposition  of  the  arc  set  into 
disjoint  subsets  Aq  exactly  as  defined  in  (2.3)  is  possible,  and 
the  definition  of  failure  pattern  X  given  in  (2.4)  also  still 
applies.  An  example  of  an  undirected  staged  network  is  given  in 
Fig.  2 . 

2.2  Recursive  Markov  Algorithms 

The  algorithms  that  we  investigate  operate  in  the  following 


way : 

(i)  A  set  of  events,  E,  has  to  be  found  where  each  event 
IeE,  is  associated  with  some  subset  of  nodes. 

(ii)  The  events  are  Markov  in  the  sense  that  their 
probabilities  of  occurrence,  pj  ,  can  be  calculated  from  the 
standard  Markov  one-step  transition  formula 


I  Pi'll, 

IeE  1  IJ 


JeE 


(2.8) 


where 


q I j  »  prob  (J  occurs)  I  occurs) 


(2.9) 


a 
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is  the  transition  probability. 

(iii)  If  the  events  can  be  numbered  so  that  transitions  go 
only  from  a  lower  to  higher  numbered  event,  then  qjj  «  0  if  J<I 
and  the  pj  may  be  calculated  recursively  from  (2.8)  in  increasing 
order  of  J.  A  special  form  of  this  is  when  the  events  can  be 

grouped  into  disjoint  sets  E1tE2 . Eq  with  events  of  one  group 

giving  rise  to  events  of  the  next  higher  group  only.  Then  (2.8) 
reduces  to 

PJ  ‘  *  PI<*IJ’  jEEq.1-  '5’1'2 . Q-  (2-'0) 

1  C  Li 

q 


3.  THE  ALGORITHMS 

3.1  A  Markov  Algorithm  for  a  Directed  Staged  Network 

Consider  a  directed  staged  network.  The  algorithm  of  Shogan 


focuses  on  events  associated  with  a  particular  stage.  Let  V  be  a 
subset  of  Sq.  Shogan  defines  the  events 


U  (V) 


I  (  V  ) 


"there  is  a  path  from  s  to  at  least  one 
node  of  V" 

"there  are  paths  from  s  to  all  nodes  of 


(3.1) 


and  gives  e xc 1  us i on- i ncl us i on  formulas  relating  the  probabilities 


of  occurrence  of  such  events  to  corresponding  events  of 
pr evi ous  stage  . 

It  is  perhaps  simpler  to  use  the  following  events. 
V^Sq_i  and  WcSq  and  define 


WWW 


[V]  =  "there  is  a  path  from  s  to  all  nodes  of  V, 

but  none  from  s  to  Sq_i\V"  (3-2) 

[V,W]  =  "Given  [V]  has  occurred,  each  node  of  W  is 
reachable  from  some  node  of  V  but  there  is 
no  path  from  Sq_i W  to  W,  or  from  V  to  Sq\W"  (3.3) 
The  event  (3.2)  can  be  viewed  as  the  union  of  a  number  of  events, 
of  a  more  general  type  given  by  Bailey  and  Kulkarni  [1],  that 
utilizes  the  staged  structure. 

Let  p[V]  and  p[V,W]  denote  the  probabilities  that  [V]  and 
[ V , W ]  occur  .  Then 

pis]  •  1 

(3.  JO 

p[w]  =  l  p[v]p[v,w],  Wc  S  ,  q  -  2 , .  .  . , Q  . 

VcS  ,  q 

-  q-  1 

This  has  form  (2.10):  the  events  {[V],  VcSq_ij  form  Eq_lf  and  the 
p [  V , Vi]  are  me  qjj.  The  last  stage  calculation  (when  q  =  Q  and 
W  =  |t})  yields  pst,  the  probability  of  s-t  connectedness. 

What  makes  (  3  •  *0  simple  is  that  the  p  [  V ,  W  ]  are  very  easily 
calculated,  because  [V,W]  depends  only  on  the  states  of  the  arcs 
of  Aq.  To  see  this,  suppose  that  veV,  weSq  and  there  is  a  path 
of  up-arcs  from  v  to  w.  From  the  definition  of  staged  network, 
this  path  may  initially  pass  through  nodes  of  Sq_i,  however,  once 
it  reaches  a  node  of  Sq  it  can  only  continue  through  nodes  of  Sq 
until  w  is  reached.  Thus,  given  [V]  has  occurred,  all  the  nodes 


of  Sq_i  on  the  path  must  belong  to  V.  If  v'  is  the  last  node  of  V 


on  the  path,  there  is,  therefore,  a  path  from  v'eV  to  w  comprised 
entirely  of  up-arcs  of  Aq.  Similarily  a  path  from  ucSq-i  W  to  W 
would,  given  that  [V]  has  occurred,  imply  the  existence  of  some 
u'eSq-iW  and  a  path  of  up-arcs  of  Aq  from  u’  to  w.  It  follows, 
therefore,  that  [V,W]  occurs  if  and  only  if,  using  paths  with 
arcs  in  Aq  only,  each  node  of  W  is  reachable  from  some  node  of  V, 
no  node  of  Sq\W  is  reachable  from  V,  no  node  of  W  is  reachable 
from  Sq_ i \ V .  This  is  the  required  result.  It  can  be  viewed  in 
the  following  way: 

L e n, m a  1  Let  V£Sq_i  and  XeFq.  Then  the  pair  V,X  Induces  a  subset 
W^Sq  defined  by 

W ( V ,  X  )  =  (w|weSq  and  there  is  a  path,  comprised  only 
of  up-arcs  of  Aq  (under  X),  from  some  node 
of  V  to  w]  (3.5) 

Q 

Figure  1  illustrates  an  example  of  such  a  W(V,X). 

Clearly 

p[V,W]  =  l  p ( X )  (3.6) 

XcF 

q 

W(V,X)  =  w 

and  ( 3  •  ^ )  reduces  to  p[s]  =  1  and 

p[W]  =  l  l  p[V]p(X),  q  -2 ,  .  . . , Q .  (3.7) 

V£S  ,  XeF 

q- 1  q 

W( V,X  )  =  W 

In  numerical  calculations  (3.7)  is  easy  to  implement  as  there  is 
no  need  to  compute  the  p[W]  one  at  a  time.  Instead,  not  3  that  in 
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calculating  all  the  p[W],  each  V,X  combination  is  used  once  only. 
Thus  (3.7)  can  be  implemented  as 
Algorithm  A 

p[s]  «-  1 
for  q  =  2  ,  .  .  ,  Q 

p[W]  «-  0  allWcSq 
for  each  V£Sq_i 

foreachXeFq  (3.8) 

W  «-  W  (  V  ,  X  ) 

p[W]  «■  p[W]  +  p[V]p(X ) 
next  X 
next  V 
next  q 

3.2  A  Markov  Algorithm  for  Undirected  Staged  Networks 

For  an  undirected  network  the  event  definition  (3*2)  does 
not  give  a  useful  algorithm  because  a  path  may  'backtrack.* 
Instead,  consider  all  the  nodes  of  a  stage  Sq  and  think  of  them 
as  parti t i oned  into  disjoint  groups.  Fig.  2  illustrates  S 2  parti¬ 
tioned  as  (1)(2)(3i4).  We  denote  a  partition  by  it  and  the  set  of 
all  partitions  of  Sq  by  Tq,  and  define  the  following  event: 

[I]  =  [it;  v] 

q 

=  " tt e T q  and,  using  arcs  of  u  A  only: 

r  =  2 

(i)  each  group  of  nodes  in  tt  is  connected 


(ii)  nodes  between  groups  are  unconnected 

(iii)  s  is  connected  to  group  v" 


(3.9) 


The  set  of  all  such  events  associated  with  Sq  will  be 
denoted  by  Eq.  Fig.  2  illustrates  the  event  [ ( 1  ) ( 2  )  ( 3 4 ) ; 3 ] e E 2 
and  the  event  [ ( 5 ) ( 6 7 ) ; 2  ]c E 3 . 

The  above  definition  can  be  viewed  as  a  specialization  of  a 
general  type  of  event  considered  by  Rosenthal  [6],  who  outlines 
an  algorithm  framework  based  on  combining  subnetworks  two  at  a 
time.  The  framework  does  not  make  specific  use  of  staged 
structure  but  clearly  has  such  a  form  in  mind.  The  algorithm 
below  takes  specific  note  of  the  Markov  nature  of  the 
calculation,  so  that  the  computational  form  is  somewhat  different 
to  that  of  Rosenthal’s  even  though  the  event  probabilities  are 
essentially  the  same.  We  use  an  analogue  of  Lemma  1: 

Lemma  2  An  event  I  =  [ir;AjeEq_i  and  failure  pattern  XeFq  induces 
a  unique  event  J  =  [p;v]eEq,  which  will  be  denoted  by  J(I,X). 

Proof  We  form  a  subnetwork,  H,  as  follows.  The  node  set  is 
Sq-iuSq.  The  nodes  of  Sq  are  treated  normally.  However,  those 
of  Sq_ 1  are  divided  by  the  partition  u  into  disjoint  groups  and, 
in  H,  the  nodes  of  each  group  are  treated  as  being  combined  into 
a  single  node  (there  being  a  separate  node  for  each  set).  The 
arc  set  of  H  is  just  the  set  of  up-arcs  of  X. 

Now  in  H  the  nodes  of  Sq  can  be  naturally  partitioned  into 
disjoint  sets;  all  the  nodes  within  a  set  are  connected  (possibly 
via  paths  which  pass  through  combined  nodes  of  Sq_-|),  but  nodes 
belonging  to  different  sets  are  disconnected.  If  we  denote  this 


partition  by  p  and  take  v  to  be  the  set  (of  p)  which  is  connected 
to  A,  then  this  uniquely  defines  J  *  [p;v].  a 

Fig.  2  gives  an  example  of  such  an  event  J(IfX). 

The  same  recursive  algorithm  ( 3  -  7 )  still  applies. 

PCs ; 1 3  =  1 

p[j]  =  l  l  p[l]p(x),  JeE  ,  q  =  2 ,  .  .  .  ,  Q 

IeE  ,  XeF  q 

q- 1  q 

JCI.X)  =  J  (3.10) 

and  again  the  modified  form  (3.8)  can  be  used  for  numerical 
calculation: 

Algorithm  B 

PCs; 1  ]  =  1 
for  q  =  2 ,  .  .  .  ,  Q 

p[J]  *  0  all  JeEq 

f or  each  I eEq_ i  (3.11) 

for  each  XeFq 

J  «■  J  ( I  ,  X  ) 

p[w]  =  p[W]  +  p[ I ]p( X  ) 
next  X 
next  I 
next  q 

4.  COMPUTATIONAL  ASPECTS  AND  COMPLEXITY 

For  Algorithm  A,  the  calculation  of  W(V,X)  and  p[W]  is  done 


times  at  stage  q.  For  Algorithm  B,  finding  the  partitions  p  is 


the  most  time  consuming  calculation.  From  the  definition  of 
J  ( I ,  X )  -  [p;v]  where  I  -  [  ir ;  A  ]  we  see  that  p  depends  on  n  and  X 
only  and  not  on  v.  Thus  in  (3-11)  p  needs  only  be  determined  for 
each  (  it  t  X )  and  not  each  (u,X,v),  a  total  of 


g( ISq-t  |  ) 


(4.2) 


times,  where  g(m)  =  number  of  possible  partitions  of  m  objects, 
is  the  sum  of  Stirling  numbers  of  the  second  kind  (see,  eg.  [5]). 
As  g  grows  hyperexponentially,  see  Table  1  for  selected  values, 
this  is  the  main  factor  limiting  the  size  of  networks  that  can  be 
handl ed . 

Table  1 


g  (m  )  = 

number 

of  partitions 

of  m 

obj  ects 

m 

1 

2 

3 

4 

5 

6 

7 

8 

2m 

2 

4 

8 

16 

32 

64 

128 

256 

g(m) 

1 

2 

5 

15 

52 

203 

877 

41  40 

Each 

calculation 

of  W 

in  Algorithm  A, 

and 

p  in  Algorithm  B, 

requires  examination  of  the  connectivity  of  the  subgraph 
Gq  =  {  Sq_ i u Sq  ,  Aq  }  .  An  algorithm  like  "breadth-first  search"  (see 
[5],  for  example)  can  be  used  for  this,  and  will  need 
0 ( | Sq_ i  |  +  | Sq | + Aq  )  elementary  operations. 

Algorithm  B  has  the  additional  requirement  of  a  hash 
function  or  subroutine  which  (i)  assigns  a  label  to  each 


partition  (of  a  set  of  given  cardinality)  and,  inversely,  (ii) 
given  a  particular  label,  specifies  the  precise  grouping  of  the 
nodes  in  that  partition.  Part  (i)  is  needed  to  identify  the 
component  of  the  array  p[J]  to  be  incremented  in  (3.11)  and  (ii) 
is  needed  to  allow  partitions  to  be  stepped  through 
systematically  in  the  outer  loop  of  (3-11).  The  Appendix  gives 
such  a  subroutine,  whose  complexity  is  0(m)  where  m  is  the 
cardinality  of  the  set  being  partitioned. 

To  test  its  effectiveness,  a  Fortran  version  of  Algorithm  B 
was  used,  on  an  IBM  Personal  Computer  AT,  to  calculated  pst  for 
the  dodecahedron  network  of  Fig.  3  and  the  grid  network  of  Fig.  4. 
The  computation  time  was  4*/2  minutes  and  52  minutes  respectively. 
Though  we  have  not  made  a  direct  comparison,  this  compares  well 
with  other  algorithms.  For  example,  for  the  test  network  used  in 
[1],  which  is  in  effect  a  dodecahedron  reduced  by  3  nodes  and  5 
arcs,  Bailey  and  Kulkarni  report  timings  of  54  minutes,  8  minutes 
and  1  minute  26  seconds  using  respectively,  the  algorithms  of 
Buzacott  [3],  Provan  and  Ball  [4]  and  their  own  [1],  on  an  IBM 
4381-k  which  is  approximately  10  times  faster  than  the  IBM  PC-AT, 
(all  algorithms  computed  the  reliability  as  .99806  when  pa  =  0.9, 
all  a )  . 

Both  Algorithms  A  and  B  resemble  dynamic  programming 
recursions  whose  calculations  of  one  stage  depend  only  on  the 
results  of  the  immediately  preceding  stage.  For  networks  like  the 
grid  of  Fig.  4,  this  means  that  a  substantial  computational 


saving  is  possible  by  avoiding  the  repetition  of  identical  calcu¬ 
lations  made  at  each  stage.  More  precisely,  suppose  the  subnet¬ 
works  Gq  -  {Sq-i uSq,Aq}  and  GP  =  { Sr_ \ u Sr  ,  Ar  }  (with  q<r,  say) 
have  exactly  the  same  structure.  Then,  during  the  calculations  at 
stage  q,  the  computed  values  of  W(V,X)  (respectively  J(I,X)) 
should  be  saved  for  all  V  cSq. j  (respectively  IeEq_i>  and  XeFq. 
This  time  consuming  computation  need  not  then  be  repeated  during 
the  rth  stage  calculation,  as,  with  suitable  labelling  of  the 
nodes  and  arcs  of  Gr,  so  that  they  correspond  to  those  of  Gq,  the 
same  values  of  W(V,X)  (respectively  J(I,X))  can  be  used  for 
V£Sr_i  (respectively  IeEr_j)  and  XeFr. 

For  example  in  the  grid  network  of  Fig.  A ,  the  CPU  time  of 
52  minutes  is  mostly  taken  up  by  the  calculations  of  stage  3  (12 
min  5  secs)  and  of  stages  ,  5  and  6  (13  min  10  secs  each).  As 
Gq  (  q  =  3  ,  ** ,  5 , 6  )  are  identical,  we  would  estimate  an  algorithm  that 
avoids  the  repeated  calculation  of  the  J(I,X)  would  reduce  the 
computation  by  more  than  half. 

Finally  we  note  that,  for  clarity  of  exposition,  an  unneces¬ 
sarily  rigid  definition  of  staging  has  been  given.  In  the  next 
section  we  show  how  the  definition  can  be  generalized.  It  should 
be  emphasized  that  though  the  timings  above  already  compare  favor¬ 
ably  with  published  work,  the  use  of  generalized  stages  can  lead 
to  further  substantial  savings  in  computing  time. 


5. 


OVERLAPPING  STAGES 


The  definitions  of  staging,  S  and  T,  can  be  relaxed  to  give 
more  flexibility  to  the  two  algorithms.  The  following  generali¬ 
zation  can  be  used  for  both  directed  and  undirected  networks. 

0:  There  are  three  sequences  of  subsets  of  N: 

S]  «  {  s } .  S2, . . . , Sq  *  { t  } 

M i  ,  M2 . Mq_!  (5.1  ) 

N  i  »  {  3  }  »  N2 . Nq  -  {t} 

such  that 

Q 

U  S  -  N  (5.2) 

q-  1  q 


Sq  “  ( Sq_ i \Mq_ i )uNq 

Mq-1  S  sq-1 
q-1 

N  n  <  U  S  )  -  <p. 
q  r-1  r 


q  =  2 ,  .  .  .  ,Q 


(5.3) 


And  for  any  ( i  , j  ) eA : 

if  ieNq,  jeNr  and  q>r  then  jeSq_i  (5. 4) 

□ 

The  condition  (5.3)  allows  the  sets  Sq  to  overlap.  Conditions 
S  and  T  comprise  the  special  case  Mq  -  Sq,  all  q,  when  Nq  -  Sq; 
the  Sq  are  then  disjoint  and  ( 5 . 4)  reduces  to  (2.2)  or  (2.6). 

It  is  readily  verified  that  the  Nq  are  disjoint  and  that 
Q 

U  N  -  N.  (5.5) 

q-1 
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The  node  stages  induce  a  staged  structure  on  the  arcs  which 

still  decompose  into  disjoint  subsets 

Q 

Ag»  Aj#...,Aq  with  U  -  A  (5.6) 

where 

Aq  -  { ( i ,j )|ieS  uNq  and  jeNq).  (5.7) 

With  this  definition  of  Sq  and  Aq,  Algorithms  A  and  B  still 
apply  completely  unaltered. 

The  stage  q  calculations  are  still  of  the  order  of  (4.1)  or 
(4. 2)  under  condition  0.  The  calculations  of  stage  q  will  thus  be 
efficient  if  Sq-i  and  Aq  are  such  that  the  following  two 
conditions  are  met. 

(i)  The  cardinality  of  Sq- j  is  kept  small  to  keep  small  the 
factor  involving  j  S  q_  ^ |  in  ( A . 1  )  or  (4.2).  This  requires  making 
| M  q_ i  j  as  large,  and  | N  q  J  as  small  as  possible. 

(ii)  The  cardinality  of  Aq  is  kept  small  to  keep  small  the 

l\l 

factor  2  in  ( A . 1  )  or  (4.2).  This  can  be  done,  indirectly,  by 

making  | N q  j  as  small  as  possible. 

There  is  a  trade-off  here  in  that  a  small  |Nq|  at  each  stage 
will  make  for  a  larger  total  number  of  stages  Q.  As  computational 
complexity  depends  only  linearly  on  Q,  but  exponentially  on  |  Aq  |  , 
the  trade  in  the  direction  of  small  |  N  q |  will  almost  always  be 
worth  making. 

As  an  example  consider  the  grid  of  Fig.  4.  If,  instead  of 
the  staging  given  in  the  figure,  we  use 


-16- 


s2  * 

(1.2 . 6) 

» 

s3  = 

(S2\(1))u7 

-  (2,3... 

,  .  ,7) 

S4  = 

(  S 3  \  (  2 )  )u  8 

=  (3.4,., 

.  .  .8) 

s5  - 

(Si|\(3)  )u9 

=  (4,5,.. 

.  .  ,9) 

30  - 

(25,26 . 

30) 

31  ' 

(t) 

then  |  S  q  |  -  6,  j  M  q_ , |  -  | Nq  J  -  1  for  q  -  2,3 . 30,  and  |Aq|  £  2 

for  q  =  3,  il,  ...30.  The  four  main  stages  of  the  original 
decomposition,  each  with  | Aq |  «  11,  is  thus  replaced  by  24  stages 
each  with  | A q |  £  2.  The  computational  complexity  is  thus  reduced 
by  a  factor  -  ( 4 x 2 1 1  ) / ( 2 4  *  2 2 )  =  85.3.  Because  of  incidental 
computational  overheads  which  we  have  not  fully  accounted  for 
(such  as  initialization  of  arrays),  the  actual  improvement  in 
speed  will  be  by  a  factor  rather  less  85.3;  however  the 
improvement  is  still  very  substantial.  A  Fortran  implementation 
of  this  version  of  staging  reduced  the  computing  time,  from  the 
52  minutes  quoted  previously,  to  1  minute  58  seconds. 

A  similar  use  of  overlapping  stages  for  the  dodecahedron 
example  of  Fig.  3  reduced  computation  time  from  4V2  minutes  to  34 
seconds  . 
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APPENDI X 

We  describe  here  a  method  of  labelling  partitions  and  its 
inverse  which  enables  a  partition  to  be  recovered  from  its  label. 

The  partitions  of  M  objects  can  be  divided  into  M  sets,  the 
kth  set  comprising  all  those  partitions  where  the  objects  are 
divided  into  k  groups.  For  instance  ( 1  3 )  (  2 )  ( ** )  represents  the 
partition  of  4  objects  into  3  groups:  (13),  (2)  and  (4);  this 

partition  thus  belongs  to  the  3rd  set.  Let 

g(k,M)  -  number  of  partitions  (of  M  objects) 
in  the  kth  set 

These  are  the  Stirling  numbers  of  the  second  kind  (see,  for 
example,  [5])  . 

We  shall  label  partitions  in  the  order  of  these  sets  (i.e. 
those  of  the  first  set  first,  then  those  of  the  second  set,  and 
so  on).  The  position  of  a  partition  is  thus  completely  defined 
once  the  position  within  the  set  to  which  it  belongs  is  fixed. 
Table  1A  illustrates  one  method  by  which  the  partitions  of  m 
objects  may  be  built  up  by  systematically  listing  all  the  ways 
that  the  mth  object  can  be  adjointed  to  each  of  the  partitions  of 
the  previous  m- 1  objects.  The  method  used  is  as  follows. 

For  a  partition  of  m  objects,  let  km  be  the  set  to  which  it 
belongs  and  jm  its  position  in  the  set.  Now,  a  partition  is 
defined  by  specifying,  for  each  object  m,  the  group  nm  to  which 
it  belongs.  In  Table  1A,  kM  and  j^  are  determined  from  nm 


(m=1,2,...,M)  using: 
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•  Hi  :  V-  v  » 


k  i  •  1  t  j  i  *  1  . 

km  “  km- 1  *  •)  m  “  ( J m- 1 “ ^) km- 1 +nm  nm  ^  km- 1 

km  -  nm,  jm  =  jm-i +g(km,m)-g(km_1 ,m-1 )  otherwise 


m»2 . M 


which  recursively  gives  the  position  jm,km  of  the  partition  of 
the  first  m  objects,  for 

Table  1A 

Partition  of  m  objects  into  k  groups 


(12) 

(123) 

(123*0 

( 1 )  (2) 

(13)(2) 

(13*0(2) 

(1  3X214) 

(1 ) (23) 

(1*0(23) 

(1 )  ( 23*4 ) 

( 1 2 )  ( 3 ) 

(12*0(3) 

( 1 2)  (3*0 
(123)(*0 

(1 ) (2) (3) 

(1  *0(2)  (3) 

(1  )(2*0(3) 

(1 )  (2)  (3*0 
(13)(2)(*0 
(1  )  (23)  ( **) 
(12)(3)(*0 

(1  )  (2)  (3)  (**) 

The  overall  position  of  a  partition  is 
k  -  1 

r  =  ""I  g(k,ra)  +  j  . 
k-1 

The  inverse  calculation  allows  the  group,  nm ,  to  which  each 
object  m  belongs,  to  be  calculated  from  r.  It  is  as  follows. 

Set  m  =  M  and  determine  km,jm  (*kM,jn)  from  the  conditions 

kM_1  kM 

r  -  I  g ( k , m )  -  jM  >  0,  r  -  £  g(k,M)  <  0. 

k  =  1  n  k-1 

Moreover  g(k,m)  =  kg(k,m)  +  g(k-l.m-l),  so  that 

A.  Ei  ther  j  m“  km8  (  km  .  m- 1  )  =  j>0,  in  which  case  the  partition  (of 
the  first  m  objects)  is  the  jth  amongst  the  last  g ( km- 1 , m- 1 ) 
partitions  of  the  kmth  set.  This  means  the  object  is  in  the  last 
group  (on  its  own!)  and  so 

nm  3  km *  km-1  =  km“ ^ »  ^m-1  *  J* 

B.  Otherwise  we  can  set  km_ j  =  km  and  find  nm,jm_1  from  the 
conditions 

nm  =  Jm  “(Jm-1-^  km  >  Jm  Jm-1km  ^  0* 

This  fixes  object  m  as  being  in  group  nm;  and  the  partition  of 

the  first  m- 1  objects  as  being  the  jm_ith  in  the  km_ith  set. 

Thus  setting  m  «•  m- 1  ,  we  can  repeat  steps  A  or  B  until  m-1, 
when  we  set  ni  -  1;  at  which  stage  all  the  nm  (m- 1  , 2 ,  .  .  .  ,  M )  have 


been  found . 
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Figure  1  A  Directed  Staged  Network 


■» 


I 

up-arc 


down-arc 


Sj  =  (s)  S2  ,1, 2, 3)  S3  =  (4,5,6)  S4  =  (t) 

A2  ~  (s,b, i • • ,g)  —  (h,i, . . . ,n)  A^  a  (o*P»r) 

Event  [V]  =  [(2,3)]  has  occurred  at  stage  2. 

A2  has  failure  pattern  X  =  (0101000). 


W(V,X)  =  (4,5). 


__________  down-arc 

Sj  »  (s)  S2  =  (1,2, 3, 4)  S3  =  (5,6,7)  S4  =  (t) 

A2  »  (a,b,...,g)  A3  =  (h,i,...,o)  A^  =  (p,q,r) 

Event  I  ■  [ ( 1) (2) (34) ; 3 ]  has  occurred  at  stage  2. 

A2  has  failure  pattern  X  =  (11010100) 


security  CLASSIFICATION  Of  THIS  RACE 


Tamm 


U  RfcfrOAT  SECURITY  CLASSIFICATION 

llnr  1  aeci  f  i  aA 


3#.  ikCuAlTV  CLASSIFICATION  AuThO«»Ty 
MA 


3b.  OCCLASSlFlCATlON/OOWNGftAOtNG  SChEOulE 

VA 


4  »>tHICHMlNG  ORGANIZATION  MtPQHt  NUMULHlSi 


6a  name  of  performing  Organization 

L'ni  versi  tv  of  North  Carolina 


REPORT  DOCUMENTATION  PAGE 


IB.  RESTRICTIVE  MARKINGS 


3.  OiSTRIBUTIONIAVAHABILITV  Of  REPORT 

Unlimited  distribution 


S.  MONITORING  ORGANISATION  REPORT  NOMUEHtSl 

AFOSR-TR-  8  6-0436 


4c.  AOORESS  (Cup.  Sl«M  and  ZIP  CmUl 

Tb.  AOORESS  (Cl ip.  SUM  and  ZIP  CixMl 

Chapel  Hill,  NC  27514 

Bldg.  410 

Bolling  AFB,  DC  20332-6448 

•a.  NAME  Of  f  UNOiNG/SPONSORING 
ORGANIZATION 


9.  PROCUREMENT  INSTRUMENT  lOENTlf  (CATION  NUMBER 


yvfosp. 


to.  SOURCE  Of  f  UNOINC  NOS. 


fROORAM 
ELEMENT  no. 


«c.  AOORESS  (Clip.  SUM  ZIP  Cod*  I 

Bldg.  410 
Bolling  AFB,  DC 


11.  TITLE  llntlmdt  security  ChaaiflcaluMI 

Markov  algorithms  for  computing  the  reliability  of  staged  networks 


13.  PtHSONAL  AUThONiSI 


WORK  UNIT 
NO. 


13a  type  Of  REPORT 

technical 


ia.  Supplementary  notation 


13b.  Time  covered 

PROM _  TO 


14.  DATE  Of  REPORT  tYt..  Ma..  Day/  IS.  PAGE  COUNT 

_ April,  1986  25 


IT  COSATI  COOES 


fifLO  1  GROUP  |  SUBGR. 


Mk  SUtijECT  TERMS  (CmImm  •«  NMfM  i^MCtMOf  Mffitifr  ly  frlocfr 

/ 

Network  reliability,  source-to-sink  connectedness, 
node  partition,  staged  network,  recursive  algorithms 


IF.  abstract  /Cunnnw*  on  m«ri<  1/  mciN«ry  mM  idfim/y  fry  frlocfr  aamfrin 

Certain  commonly  occurring  types  of  network,  whether  directed  or  undirected, 
exhibit  a  staged  structure.  Two  algorithms,  based  on  node  partitioning, 
are  presented  which  take  advantage  of  such  structure  and  which  use  a  Markov 
transition  probability  form  of  recursion.  The  algorithm  for  directed  networks 
is  related  to  the  Markov  chain  formulation  of  Bailey  and  Kulkarni,  but 
for  undirected  networks  a  more  detailed  form  of  state  definition  is  used 
related  to  one  suggested  by  Rosenthal. 

1 

1  The  computational  advantages  of  the  algorithms  are  discussed  and  some 

J  numerical  results  presented. 


30.  OlETRiauTlON/AVAILAEILlTV  OP  ABSTRACT 
UNCLASSlfiCO/UNLIMlTSO  Q  SAME  AE  RPT.  □  OTIC  USSR*  □ 


33*.  NAME  OP  RESPONSIBLE  INOlVIOUAL 


31.  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


S3*  OffiCC  SYMBOL 


lit  TELEPHONE  NUMBER 
HtlltM  Aura  CM«| 


_  R.C.H.  Chen 


00  FORM  1473.  83  APR 


niiWFi;i. 


EDITION  OP  I  IAN  U  IS  OBSOLETE. 


security  Classification  op  TpiiS  pace 


